A young protoplanetary disk in the Bok globule CB 26? 
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ABSTRACT 

We present sub-arcsecond resolution millimeter-wave images of a circumstellar disk in the Bok 
globule CB 26. The presence of an edge-on disk is confirmed by the dust continuum morphology 
and the velocity field of ^^CO emission, which displays a Keplerian rotation pattern about an axis 
perpendicular to the long axis of the dust emission. We deduce a mass ^O.SM© for the obscured 
central star. The disk is optically thick at mm wavelengths inside 120 AU, has a symmetric 20° 
warp beyond 120 AU, an outer radius of ~ 200 AU, and a mass of at least 0.1 Mq. We suggest 
that the CB 26 system is in an intermediate stage between deeply embedded protostellar accretion 
disks and the more evolved, perhaps protoplanetary, disks around T Tauri stars. 



Subject headings: circumstellar matter 
toplanetary disks — stars: formation 

Introduction 



dust, extinction — ISM: globules — planetary systems: pro- 



CB 26 (L 1439) is a small cometary-shaped Bok 
globule located '^10° north of the Taurus/ Auriga 
dark cloud. Launhardt & Henning (1997) sug- 
gested a distance of 300 pc. From a re-examination 
of the larger-scale velocity structure (Ungerechts 
& Thaddeus 1987), we believe it is part of the Tau- 
rus/Auriga complex at 140 pc. Single-dish submm 
images at 10" resolution show strong imresolved 
dust continuum emission at the south-west rim of 
the globule, surrounded by a thin asymmetric en- 
velope with a well-ordered magnetic field directed 
along PA'- 25° (Henning et al. 2001). A small 
bipolar near-infrared (NIR) nebula bisected by a 
dark extinction lane is associated with the millime- 
ter emission. The sub-arcsecond resolution NIR 
polarization pattern is consistent with the pres- 
ence of an almost edge-on circumstellar disk elon- 
gated along the extinction lane at PA ^ 60° and a 
thin envelope (Stecklum et al. 2001). The spectral 
energy distribution indicates a low-mass, > 0.5 L0 
pre-main sequence (PMS) star surrounded by a 
disk and an envelope (Stecklum et al. 2001), i.e., a 
Class I object (Adams, Lada, & Shu 1987; Kenyon 
& Hartmann 1995). 



2. Observations and data reduction 

CB 26 was observed with the Owens Valley Ra- 
dio Observatory (OVRO) millimeter-wave array 
between January and December 2000. Three 
configurations of the six 10.4 m antennas pro- 
vided baselines in the range 6-180kA at 2.7 mm 
(llOGHz) and 12-400kA at 1.3mm (232GHz). 
Average SSB system temperatures of the SIS re- 
ceivers were 300- 400 K at 110 GHz and 300- 
600 K at 236 GHz. The digital correlator was 
centered on the "CO(l-O) line at 110.2 GHz, 
adopting the systemic velocity of CB26, ulsr = 
5.5kms^^; spectral resolution and bandwidth 
were 0.17kms~^ and 5kms~^, respectively. Con- 
tinuum emission at 1.3 and 2.7 mm was observed 
simultaneously in 2 GHz-wide bands, except in the 
highest-resolution configuration, where the 4 GHz 
wide-band capability of new 1 mm receivers was 
used. Amplitude and phase calibration relied on 
frequent observations of a nearby quasar, result- 
ing in absolute position uncertainty of 0'.'2. Flux 
densities are based on observations of Uranus and 
Neptune, with resulting uncertainties of 20%. The 
raw data were calibrated and edited using the 
MMA software package (Scoville et al. 1993). 
Mapping and data analysis used the MIRIAD tool- 
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box (Sault, Thcubcn, & Wright 1995). Wc also 
applied maximum entropy deconvolution (MEM) 
to the data and obtained images which are very 
similar to those obtained by standard cleaning. 

3. Results 

In Fig. 1, our dust continuum images of CB 26 
show a source elongated along PA = 60 ±5°. At 
2.7 mm, the projected length derived from the dis- 
tribution of clean components is 230±30AU. The 
minor axis is unresolved and there is no hint of an 
envelope. By contrast, the 1.3 mm images show 
a narrow central structure of length 400±40AU 
and a small envelope of size 220 AU x 280 AU. 
The scale height h of the central structure re- 
mains unresolved, i.e., /i<20AU. Beyond radii 
of 120 AU, an almost symmetric 20° warp is ap- 
parent in Fig. Id. Total continuum flux den- 
sities are 22±4mJy and 150±30mJy at 2.7 and 
1.3 mm, respectively. Thus, we recover completely 
the 1.3 mm flux of the unresolved component as 
measured with the IRAM 30-m telescope, 160 mJy 
(Launhardt & Henning 1997). From our maps, we 
calculate that about 50 mJy arises in the envelope 
and 100 mJy in the narrow central structure. Of 
this 100 mJy, about 75mJy derives from the inner 
230 AU where the 3 mm emission arises. 

In Fig. 2, the 1.3 mm continuum MEM image 
is overlaid on a map of the NIR reflection neb- 
ula. Within the NIR astrometric uncertainties, 
0'.'5, the dust emission coincides with the extinc- 
tion lane and with the location and orientation 
of the disk postulated by Stecklum et al. (2001). 
Our kinematic data support the disk hypothesis. 
As illustrated in Fig. 3a, strong ^^CO(l-O) emis- 
sion (T^ 20 K) arises mainly from the outer 
parts of the disk, suggesting optically thick emis- 
sion. In Fig. 3b, the ^^CO line is clearly double- 
peaked, with blue-shifted emission concentrated 
to the north-east and red-shifted to the south- 
west. The velocity structure indicates a rotation 
axis perpendicular to the projected plane of the 
disk. At the systemic velocity, ulsr = 5.5 kms^^, 
gas emission appears to be self-absorbed, probably 
due to the extended outer envelope. 

4. Discussion 

Our millimeter interferometric observations of 
thermal dust continuum and ^^CO line emission 



a) 2.7mm 
^^^1 1 1 1 1 1 1 ■ ' ' ^ t 1 U 1 1 1 1 ^ 


c) 2.7mm (MEM) 

■ % ; 

' ' 1 1 1 1 1 1 1 ' 1 ^ ' 1 ' 1 1 1 1 1 1 M 1 1 


b) 1.3mm ' - 


d) 1.3mm (MEM) 















2 1 0-1-2 2 1 0-1-2 
AR.A. ["] AR.A. ["] 



Fig. 1. — Dust continuum emission from CB26 
(R.A.=04:59:50.74, DEC=-h52:04:43.8, J2000): 
la) and b) Cleaned and restored 2.7 and 1.3mm 
maps. Contour levels arc at -2, 2, 4, 7 ... 
22 X Icrrms where a = 0.7mJy/beam at 2.7mm 
and 1.3mJy/beam at 1.3 mm. Beam sizes are 
shown as grey ovals (l'.'2x0'.'84 at 2.7mm and 
0'.'58x0'.'39 at 1.3mm). c) and d) Deconvolved 
maps derived using the maximum entropy algo- 
rithm (MEM). 
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Fig. 2. — MEM image of the 1.3 mm dust emis- 
sion (solid white contours at 4, 11, 18, 25, and 
32mJy/arcsec^) overlaid on a J — H color map of 
the bipolar NIR reflection nebula. Black contours 
show the nebula K-band emission. The white cross 
denotes the presumed location of the illuminating 
star (from Stecklum et al. 2001). 
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Fig. 3. — a) Integrated ^"^00(1 0) emission from 
CB26 in grey-scale (dotted contours at 8, 13.5, 
19Kkms~^) with 1.3mm dust continuum con- 
tours overlaid. A dashed line represents the plane 
of the disk. The CO beam size is shown as grey 
ellipse, b) Position-velocity diagram along the 
disk major axis. Thick contours at 27, 45, 63, 81, 
99% of the maximum intensity show the observed 
velocity field. The expected emission from a Ke- 
plerian disk around a 0.35 M0 star is represented 
by dashed contours (see text). 

from CB 26 confirm the presence of an edge-on 
circumstellar disk as postulated by Stecklum et al. 
(2001). The disk has a diameter of about 400 AU 
and an as yet unexplained warp beyond radii of 
120 AU. This could be due to an undetected wide 
companion star with an inclined orbit. 

The continuum fluxes lead to an average 1.3- 
3 mm dust emission spectral index a for the in- 
ner i?< 115 AU of the disk of 1.7±0.2, whore flux 
Sn (X v'^ . Since a = 2 for a blackbody in the 
Ray leigh- Jeans limit, emission from the inner disk 
must be very optically thick, making it difficult 
to constrain the mass. Beyond 120 AU, no 3 mm 
emission is detectable and a > 2.7±0.3, suggest- 
ing optically thin emission. The transition be- 
tween optically thick and optically thin emission 
occurs approximately at the radius where a warp 
appears, ~120AU. Thus the outer disk may be 
disturbed and inflated, and flared T Tauri star 



(TTS) disk models may no longer apply. Nev- 
ertheless we derived a lower limit to the disk 
mass by comparing the 1.3 and 2.7mm contin- 
uum emission to spatially resolved radiative trans- 
fer calculations of different heuristic 'typical' TTS 
disk models characterized by temperature profiles, 
Td oc r"*^, with 5 = 0.5-0.75, surface density pro- 
files, S oc r~P, with j3 = 1.0-1.5, and disk flar- 
ing, H oc r^, with 2; = 0-1.5 (e.g., Adams et al. 
1987; Kenyon & Hartmann 1987). An interstel- 
lar gas-to-dust mass ratio of 100 and typical disk 
dust opacity Ki.3 = 2cm^g~^ with k oc and 
(3 — 1 were assumed (cf. Beckwith et al. 1990). 
All models suggest disk masses of at least 0.1 Mq. 
Due to the high optical depths of the ^^CO line, 
self-absorption, and possible depletion (cf. Thi et 
al. 2001), an independent mass estimate from the 
molecular line emission is unobtainable. 

The envelope detected in the 1.3 mm images 
(Figs, lb, Id) is probably associated with an 
outflow or disk wind. Assuming optically thin 
emission at temperature 30 K and K1.3 = 1 ± 
0.5cm^g~^, a fairly typical value for dense proto- 
stellar cores (cf. Ossenkopf & Henning 1994), we 
derive MH(env) = 0.01±0.003Mq. For a cylindri- 
cal morphology (cf. Figs, lb. Id, and Sect. 3), we 
infer an average density (uh) ~ 310^cm~^, col- 
umn density Nn ~ 510^'^cm~^, and correspond- 
ing visual extinction Ay ~ 300 mag, implying that 
back- warming from the envelope has to be consid- 
ered in the energy balance of the disk (cf. But- 
ner, Natta, & Evans 1994). The total mass of the 
more extended (~ 3000 AU diameter) asymmetric 
envelope seen in single-dish maps is 0.1 ± 0.05 Mq 
(Henning et al. 2001). This may be a remnant of 
the globule core from which the system formed. 

The velocity structure of ^"'CO can be modeled 
by Keplerian rotation when the high optical depth 
of the inner disk and self-absorption due to an ex- 
tended envelope are taken into account, yielding a 
mass of the central star of O.35±O.1M0 (cf. Fig. 
3b). However, since the disk mass (>O.1M0) is 
comparable to the inferred mass of the central star, 
self-gravitation may be important in the outer disk 
(cf. Bertin & Lodato 1999). The outer rotation 
curve would then be flatter than in the pure Keple- 
rian case, and the mass of the central star concomi- 
tantly lower, perhaps only 0.25 Mq. The ^^CO line 
strength suggests a gas-rich disk with kinetic tem- 
perature >20K outside R ~ 100 AU. An upper 
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limit to the turbulent velocity dispersion in the 
disk is O.Skms"^. 

PMS evolutionary tracks indicate an age of 
7 ± 2 X lO^yr for a 0.3 Mq star with L = 0.5 
(D'Antona & MazziteUi 1994). Due to the 
anisotropic radiation field caused by the optically 
thick edge-on disk, 0.5 Lq may be only a lower 
limit to the bolometric luminosity of the central 
star (see Men'shchikov & Henning 1999). But, for 
such a massive disk, the intrinsic viscous luminos- 
ity could also contribute. Overall, the properties 
of the CB 26 system are consistent with its being 

younger than classical TTSs (ages 10® 10 ''yr). 

There is no dense, centrally peaked cloud core with 
evidence of collapse, and no prominent molecular 
outflow, indicating that the main accretion phase, 
assumed to last a few times 10** yr (e.g., Andre, 
Ward- Thompson, & Barsony 2000), has ended. 
However, the disk mass is close to the theoreti- 
cal gravitationally unstable limit, 0.3 M*, so that 
the envelope may still be accreting at low rates 
onto the disk (see HoUenbach, Yorke, & John- 
stone 2000, and references therein). Indeed, Fig. 
3b shows some 'forbidden' red-shifted emission 
from the blue side of the disk which may be due 
to infall. 

We conclude that the CB26 disk surrounds 
a PMS M-type star with age ~10^yr. This 
represents an intermediate stage between such 
deeply embedded protostellar accretion disks as 
L1551IRS5 (e.g., Butner et al. 1994), and 
the more evolved, perhaps protoplanetary, disks 
around classical TTSs, which may be the reser- 
voirs for material for planet formation. 
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